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Abstract 

It is well known that in gauge theories, the spin (and orbital angular momentum) of gauge 
particles is not gauge invariant, although the helicity is; neither are the canonical momentum and 
canonical angular momentum of charged particles. However, the simple appeal of these concepts 
has motivated repeated attempts to resurrect them as physical descriptions of gauge systems. In 
particular, measurability of the gluon-spin-contribution to the proton helicity in polarized proton 
scattering has generated many theoretical efforts in generalizing it and others as gauge-invariant 
quantities. In this work, we analyze the constraints of gauge symmetry, the significance of gluon 
spin in the light-cone gauge, and what is possible and natural in QCD parton physics, emphasizing 
experimental observability and physical interpretation in the structure of bound states. We also 
comment on the measurability of the orbital angular momentum of the Laguerre-Gaussian laser 
modes in optics. 



1 



I. INTRODUCTION 



Feynman's parton is a powerful idea describing high-energy collisions of hadrons [l[ . The 
parton concept grows out of a simple model for hadrons in the infinite momentum frame 
(IMF), i.e., observers or systems traveling at the speed of light, in which the interactions 
are slowed down by Lorentz dilation and the hadrons appear to be simply a collection of 
non-interacting partons. When quantum chromodynamics (QCD) was proposed as the fun- 
damental quantum theory to describe the physics of strong interactions, the parton picture 
survived with two minor modifications: First, it emerges in a special choice of gauge, the 
light-cone gauge A + = (where + means a combination (A + A 3 ) / \/2 of A^, /i = 0, 1, 2, 3, 
if the hadron is moving in the 3-direction), because the quarks do interact with each other 
through an SU(3) color gauge field, even in the IMF. In this choice of gauge, the interaction 
effects are nominally gone. Second, partons act differently at different (momentum) resolu- 
tion scales, again due to color interactions. However, the scale evolution can be computed in 
QCD perturbation theory and is believed to be under control [2j. Of course, QCD also per- 
mits studying effects beyond the simple parton model by including parton correlations, trans- 
verse momentum, off-shellnesss, etc. These effects cannot be studied model-independently, 
and one must invoke high- twist factorization and expansions 0, H[ . 

Because the parton concept is so appealing, it often motivates attempts to do away with 
the gauge symmetry of QCD and search for "physical observables" , which are intuitive in 
non-interacting theory but "unnatural" due to gauge symmetry constraints. These attempts 
have occasionally met remarkable successes. For example, although the gluon spin is not a 
gauge- invariant concept, the gluons in the IMF have momentum predominantly along the 
3-direction and appear to be on-shell. A free gluon is defined to have helicity of either 
+ 1 or —1. It is then a simple matter to define the gluon helicity density Ag(x) as the 
density difference of two different helicities. The total gluon polarization AG = J Ag(x) 
is a measurable quantity through hard scattering, and therefore must be gauge invariant. 
On the other hand, in the rest frame of the nucleon, the transverse component of the gluon 
momentum is by no means negligible, and the projection of the gluon spin along the 3- 
direction is no longer gauge invariant. Moreover, as a part of the bound state, the gluons 
are definitely off-shell. An off-shell gluon, as in the gluon propagator, is hard to define 
gauge-invariantly. Therefore, Ag(x) can no longer be interpreted simply as the bound state 
property of the nucleon, although it is a matrix element of a gauge-invariant operator, 
nonlocal and complicated in general framework [Bj]. Nevertheless, in the light-cone gauge, 
it reduces to the simple gluon spin operator along the 3-direction, reminiscent of the IMF 
physics. 

The example of the IMF gluon helicity has motivated myriad attempts to define parton 
momentum and orbital angular momentum (OAM). In particular, there has been a recent 
attempt to completely reinvent the concept of gauge invariance by proposing to decompose 
the gauge potential into the so-called physical and pure-gauge parts, and defining various 
"physical observables" with the physical part 0,0|. The standard textbook result that gluons 
carry about one-half of the nucleon momentum has been challenged by the new definition 
of the parton canonical momentum. There have been a number of follow-up studies in the 
literature [8-12 



In a different context, there have been many studies about the parton transverse- 
momentum- dependent (TMD) factorization "theorems" in the literature. In some cases, 
these theorems can be argued or even proved and the relevant distributions are gauge invari- 
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ant [131 ] . In many other cases, however, such factorizations are simply model assumptions, 
and do not meet the requirement of gauge symmetry. In particular, when the parton trans- 
verse momentum is fully taken into account, the twist-expansion is often lost in the sense 
that all twists are now involved. There seems to be no systematic expansion since there are 
no arguments as to why extra gluons with transverse polarizations should not be included 
in the calculations (l4j . 

In this paper, we analyze the constraints of gauge symmetry on bound state properties 
and parton physics. In Sec. II, we make some general remarks about gauge symmetry, 
particularly on how to maintain gauge symmetry in a fixed-gauge calculation. We emphasize 
that the wave function is a gauge- dependent concept and that a gauge symmetric-matrix 
element emerges only when the operator is manifestly gauge symmetric, modulo the so- 
called BRST exact operators and the equations-of-motion operators in the covariant gauge. 
In Sec. Ill, we analyze the gauge symmetry in parton physics, particularly in the case 
of gluon polarization. We point out that while gluons appear as on-shell radiation in the 
IMF, they are off-shell and have transverse momentum as a part of the bound state. We 
introduce the concept of gauge-invariant extension (GIE) to establish a connection between 
parton physics and gauge symmetry. In Sec. IV, we discuss the GIE of the Coulomb gauge 
in connection with the work of Chen et al. We comment that the result of Wakamatsu's 
calculation of the scale evolution for the gluon spin is in fact a light-cone gauge result 
We will also discuss the spin and OAM in the case of photon radiation fields in optics. We 
conclude the paper in Sec. V. 



II. GENERAL REMARKS ABOUT GAUGE SYMMETRY 

In this section, we make some general remarks about gauge symmetry. We emphasize 
that wave functions in non-relativistic quantum mechanics are gauge dependent, and the me- 
chanical momentum of charged particles corresponds to the covariant derivative in quantum 
theory. We argue that in atomic systems, because of the non-relativistic approximation, the 
partial derivative is actually gauge invariant to leading order inv/c expansion. We then con- 
sider the gauge symmetry in quantum field theories, reminding the reader why gauge choice 
is a necessity and how a gauge-invariant result emerges from a fixed-gauge calculation. We 
also comment on the relationship between gauge symmetry and Lorentz symmetry. These 
remarks are by no means novel, and many are present in textbooks or papers. However, 
they lay out an important general framework for gauge symmetry and will be useful for the 
discussions in the remainder of the paper. 



A. Gauge symmetry in non-relativistic quantum mechanics 

Let's first consider gauge symmetry in non-relativistic quantum theory with electromag- 
netic interactions. For a charged particle moving in an external electromagnetic field, the 
Hamiltonian is, 

where P is the canonical momentum and we have customarily called A = <ft. It has the 
following eigenvalue problem, 

H^ n {r) = EMr) , (2) 
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with energy eigenvalues E n and energy eigen functions ^ n {r). 
Under a time-independent gauge transformation, 

A"(r) y4 M '(f) = A"(r) + <9 M x(0 

we obtain a new Hamiltonian, 

2m 



(3) 



(4) 



which is manifestly different. However, it has the same quantum mechanical eigen energy, 

H%(f) = E n ^' n {r) , (5) 

where ip' n = e iex ^ip n (f) is the eigen function in the new gauge. Thus, while the charged 
particle probability density, 

Pn{r) = <(r)^n(0 , ( 6 ) 

is gauge invariant, the expectation value of the charged particle canonical momentum P 
(which is quantized clS db partial derivative) is not: 



tf m P = U m P ^ ) + e |V*(r)l ^„> . 



(7) 



Thus, the proof that the canonical momentum matrix element is gauge invariant in Ref. 



11] is likely incorrect. Actually, in Ref. the gauge transformation leaves the physical 



states invariant, but this is not true according to Eq. (5): wave function is a not a gauge- 
independent quantity. It is the matrix element of the mechanical momentum p = P — eA 
that is gauge invariant: 



1L 



P-eA' 



P -eA 



(8) 



which corresponds to the covariant derivative in quantization, p = P — eA = —iD. 

In Feynman's lectures on physics, he gave a simple example to demonstrate that it is the 
covariant derivative that corresponds to the mechanical momentum in quantum mechan- 
ics [TBI]: Consider a charged particle near a solenoid which has a quantum mechanical wave 
function ip(r,t). The solenoid does not have any current in the beginning. At some point in 
time, a current passes through the solenoid and a stable magnetic field is established. In the 
process, the particle gets a momentum kick because a changing magnetic field induces an 
electric field which exerts a force on the charged particle. From the Schrodinger equation, 
we know that the wave function must be continuous in time. Therefore, the momentum kick 
on the particle cannot be obtained from the partial derivative acting on the wave function, 
which must be also continuous. Instead, this momentum kick comes from the establishment 
of the vector potential A in the system. 



B. Atom systems and the Coulomb gauge 

Atoms are, to a very good approximation, non-relativistic systems. In such systems, the 
charge effects are much larger than current interactions, and the dominant interaction is 
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the static Coulomb force. Let us establish a power counting in terms of the velocity v and 
electromagnetic coupling e [l6j]. The fine structure constant e 2 /4ir = small quantity, 

and so is the typical velocity of the charged particle v/c ~ a em . We will show through 
power counting that the canonical momentum and canonical OAM are approximately gauge 
invariant. 

In non-relativistic systems, 

B/E ~ v/c . (9) 

We are going to use this to establish a power counting for the four-vector potentials. Ac- 
cording to the definition: 

dA 

E = -V<j>-—, B=VxA, (10) 

since the electrostatic interaction e<p is order of a 2 m (recall the Bohr energy), we have 

<f) ~ ea cm ; E ~ eva em . (11) 

In principle, we can choose A to be of the same order of magnitude as 0. However, since B 
is suppressed by a factor of v/c relative to E, the transverse part of A must be suppressed 
by the same factor. To have a natural power counting, one should work only in the class of 
gauges in which the longitudinal part of A is similarly suppressed. One gauge that satisfies 
this condition is the Coulomb gauge, 

V-A = 0. (12) 

which eliminates the longitudinal part of the gauge potential entirely. Thus the Coulomb 
gauge is a natural choice for non-relativistic systems. 

From B/E ~ v/c, the counting for the vector potential in Coulomb gauge is 

A ~ eva ern . (13) 

Thus eA in the Coulomb gauge is a 2 m suppressed relative to the canonical momentum 
P = —iV ~ v. Therefore, 

p = P + 0(va 2 cm ) , (14) 

and the canonical momentum is gauge invariant to the relative order a 2 m . This explains 
why the canonical momentum and canonical OAM are useful concepts in atomic systems. 

Of course, one can choose to work in a gauge in which eA is not suppressed at all relative 
to P. Then the atomic wave functions will be different, and the physical momentum and 
OAM must follow from the mechanical momentum and mechanical OAM. 

In QCD bound states, such as the nucleon, the above power counting fails. In particular, 
the color magnetic field B a is not suppressed relative to the color electric field E a . Moreover, 
the quark motion is entirely relativistic. Therefore, the eA and P have the same order of 
magnitude. The Coulomb gauge is no longer special to keep the power counting natural. 
In factorization theorems involving parton transverse momentum, it is often assumed that 
parton canonical momentum is much larger than the transverse gauge potential in the light- 
cone gauge. There is no known physical basis for this assumption. 



5 



C. Quantization of QED and gauge symmetry 



Here we consider gauge- symmetry constraints in quantum field theories. For simplicity, 
we take the U(l) gauge theory with a spin- 1/2 fermion, such as quantum electrodynamics 
(QED), as an example. The Lagrangian density is 

C = ^(i^D* - rr# - \f^F^ , (15) 

where ip is the electron's Dirac field and is the U(l) gauge potential. 

has two physical degrees of freedom; it also mediates static electric and magnetic 
interactions among charged particles. However, N 1 does have spurious degrees of freedom 
that must be constrained in order to have a finite photon propagator. Thus, quantization of 
a gauge theory requires gauge fixing, which ironically breaks the gauge symmetry. Although 
the symmetry is no longer manifest after quantization, the physical observables in terms 
of the physical degrees of freedom shall be independent of the spurious gauge dynamics. 
This can be verified either through residual gauge symmetry after gauge fixing or through 
comparing calculations in different gauges. 

Different gauge fixings often lead to different Hilbert spaces, and the physical states are 
gauge dependent, just as in the case of non-relativistic systems described in the previous 
subsection. For instance, in the Coulomb gauge, all quanta in the Hilbert space are physical 
degrees of freedoms, whereas in the Weyl gauge A = 0, the longitudinal polarization of 
the photon is allowed, and the physical states are chosen through Coulomb's law constraint 
(V ■ E — p)|^ / p hys) = 0. In the covariant Lorentz gauge, all four polarizations of the pho- 
tons are allowed and the physical Hilbert space is constrained by Lorentz gauge condition 
c^A^l^phys) = 0. However, due to the residual gauge freedoms, the physical state is still not 
unique. It can contain non-physical quanta. 

Despite the fact that all Green's functions and wave functions are gauge dependent, the 
poles of the Green's functions and scattering S-matrix are gauge invariant. Usually, the 
proof is done using residual gauge symmetry. Here we are concerned with the physical 



matrix elements in asymptotic states (ip\0\(f)) (11] . In order for them to be gauge invariant 



we have to study general gauge symmetry transformation. Let's assume a superoperator G 
(in the sense that it is capable of transforming states from one Hilbert space to other states 
in a different one) which brings the physical states from one gauge to another, 

W) = G\f) . (16) 

Then the matrix element is invariant under gauge transformation, 

WPWW = Q>\0{A)\4>) , (17) 

only if O(A') = 0(A) is gauge invariant; that is, it has the same dependence on the gauge 
potential as in O. 

The requirement for matrix elements to be gauge invariant is much simpler to see in the 



path-integral formulation [17] . The sufficient condition again is that the operator has to be 
manifestly gauge invariant. 

Of course, it is not necessary to have gauge-invariant operators in order for its phys- 
ical matrix elements to be gauge invariant. In covariant quantization, it is known that 
the equations-of-motion operators and the so-called BRST-exact operator do vanish in the 
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physical states 18[. For an example of the energy-momentum tensor, see Ref 19|. It is 
also taken as true that the gauge-noninvariant total derivative operators do vanish in well- 
localized states. However, in plane wave states, the total derivative operators do contribute 
to the non-forward matrix elements. Thus one has to be careful with the total derivative 
operators, particularly when they are gauge dependent. 

The matrix element of a gauge- dependent operator usually vanishes when averaged over 
all gauges, just like the velocity vector vanishes when summed over contributions from all 
frames 20 



D. Lorentz symmetry 

According to Einstein's theory of special relativity, all physical quantities should be 
Lorentz tensors, i.e., when coordinates transform, they behave either as scalars, four- vectors, 
or high-order tensors. Therefore, when making a frame transformation, one easily obtains 
the physical quantities in one frame from another. In gauge theories, because not all four 
components of A M are physical, calculations are found not to be manifestly Lorentz covariant. 
For example, the gauge choices often break Lorentz symmetry when the gauge conditions 
are frame dependent. However, for local gauge-invariant quantities, Lorentz symmetry is 



manifestly kept in calculations [21|. For nonlocal operators, one has to introduce gauge 
links, which often involve geometrical lines. The transformation of such a line is usually 
complicated, and Lorentz symmetry transformation of nonlocal operators often becomes in- 
tractable. In this case, one considers the subclass of Lorentz transformation which leaves 
the line invariant. 

In the case of noncovariant gauges, calculations in different frames often correspond 
to different gauge choices. Therefore, gauge symmetry of the calculations is a necessary 
condition to maintain Lorentz symmetry, or vice versa. We note that in quantum theory, 
the wave function of a system is not a Lorentz invariant quantity — it changes with the 
coordinates. Such a change is not just a kinematic transformation; it involves a Lorentz 
boost which is usually interaction dependent. Thus boosting a wave function is a highly 
nontrivial matter in a relativistic interacting theory. 



III. EXTENDED GAUGE INVARIANCE IN PARTON PHYSICS 

In high-energy scattering in QCD, partons are often used to describe the scattering pro- 
cesses involving large momentum transfers Total scattering amplitudes usually contain 
two factorized parts, described as factorization theorems [1] . One part contains the product 
of parton densities and/or fragmentation functions and/or soft functions. As far as the hard 
scattering is concerned, the hadrons can be viewed as beams of non-interacting partons, 
and the parton densities simply describe the properties of the beams. The second part is 
the parton scattering cross section, which involves scattering of on-shell colored particles. 
Feynman's parton language greatly simplifies the description of high-energy scattering. The 
proof of the factorization theorems is to show that the cross section has a systematic expan- 
sion in terms of the inverse powers of large momenta and that the leading factorized parts 
are naturally gauge invariant. 
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A. Gauge symmetry of parton density and physics in the light-cone gauge 

One of the keys to the parton language description of high-energy scattering is gauge 
symmetry and the light-cone gauge, A + = 0. [Here we use the light-cone coordinates 
A ± = (A ±A 3 )/\/2, A±_ = (A 1 , A 2 ), and it is well known that the light-cone coordinates are 
equivalent to the use of IMF.] The parton density can be described as the matrix element 
of gauge-invariant operators in hadron states, and since the cross section is gauge invariant 
by itself, this means that the on-shell parton scattering cross section is also gauge invariant 
and can be calculated independently in any gauge. It can also be seen that the parton 
cross sections involve only on-shell parton scattering. The gauge invariance, however, masks 
the physical property of the parton density operator. For example, the quark longitudinal- 
momentum density distribution, 

q(x) = ^J ^(PSmXn^mPS) , (18) 

where n is a light-like four- vector n 2 = and the nucleon momentum is P M = (P°, 0, 0, P 3 ), 
with P ■ n = 1. is a gauge- invariant quark field defined through multiplication of a 

light-cone gauge link, 

= exp \-ig £ n ■ A(i + Xn)dXj V>(£) . (19) 

This gauge link ensures that whenever a partial derivative or canonical momentum of colored 
quarks appear, the gauge potential A^ must be present simultaneously to make it a covariant 
derivative (mechanical momentum), = <9 M + igA^ . Indeed, taking its moments, 

J x n - 1 q(x)dx ~ n^ 1 ...n^(P[4;(0)Y 1 iD fl2 ...iD fl "tp(Q)\P) , (20) 

we see that the parton momentum distribution refers to the gauge-invariant mechanical 
momentum! The mechanical momentum structure is clearly seen through Feynman diagrams 
in Fig. [T] Gauge symmetry requires that a parton with mechanical momentum k + = xP + 
includes the sum of all diagrams with towers of longitudinal gluon A + insertions. 




FIG. 1: Deep-inelastic scattering process in which the gauge invariance involving the longitudi- 
nal quark mechanical momentum xP + is achieved through insertions of gluons with longitudinal 
polarization A + . 

Simple parton physics emerges in the light-cone gauge A + = 0, where the light-cone 
gauge link disappears and all the covariant derivatives in Eq. (120]) become partial ones, 
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i.e., the canonical momentum and the simple quark field become physical. One can use the 
light-cone quantization to write: 



MO 



dk + d 2 k i 



d ] (k + , k ± )v{k + , k ± )e^ + b{k + , k ± )u{k + , k ± )e- iik+ t~~^<^ 



J 2A;+(2tt) 3 

Inserting the above expression into the parton density, one finds for x > 0, 



q(x) 



d 2 k. 



2k+(2ir) 



[P\tf(xP + ,k ± )b(xP + ,k ± )\P) 



(21) 
(22) 



which is a genuine number-density distribution. 

While the parton language simplifies the description of a scattering process considerably, 
one does not have to choose the light-cone gauge A + = 0. There always exists a natural 
gauge-invariant formalism that embodies the factorization, such as the presence of the matrix 
elements of the twist-two gauge-invariant operators. 

In recent years, there have been many discussions on the so-called /c^-dependent factor- 
izations. The motivation is to take into account the effects of parton transverse momentum 
in the scattering processes. In some cases where there is a systematic power counting for 



the relevant physical observable, the factorization is gauge invariant [13] . However, in many 



other cases, the factorization is basically a model assumption which cannot be justified be- 
cause other contributions of the same order have been neglected and the result is not gauge 
invariant 14 . 



B. Gluon spin 

It is well known that in gauge theories, the spin operator of a gauge particle, defined as 

S g = J d 3 £E x A , (23) 

is not a gauge- invariant quantity [22[, because of the manifest dependence on the gauge 
potential A. This applies for photons in QED and gluons in QCD. However, its helicity, 
i.e., the projection of S g along the direction of momentum, is gauge invariant. This can be 

easily seen from the vanishing of the gauge-dependent term (V x V%) • E. Helicity is also a 
projection of the total angular momentum along the direction of motion, as the orbital part 
obviously does not contribute. The helicity of a vector gauge particle is defined to be +1 or 
— 1, corresponding to right-handed or left-handed polarization, respectively. 

While helicity is a useful gauge- invariant concept, it should not be confused with spin 
projection. For example, the spin vector S g calculated along the _L-direction for a photon 
moving in the 3-direction definitely depends on the choice of the wave function or gauge. 
Many concepts developed in radiation physics use helicity as a starting point. For example, 
the well-known electromagnetic multipoles are constructed with Wigner rotation matrices 
applying to plane wave states with good helicities j22[ . 

A nucleon with large 3-momentum P 3 — > oo is made of off-shell (virtual) constituents 
with large 3-momenta, xP 3 , as well. In high-energy scattering, only this 3-momentum is 
important and the transverse component and off-shellness can be neglected because they 
produce subleading effects in the parton scattering cross section. Therefore, as far as the 
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hard scattering process is concerned, the nucleon appears to be made of a beam of massless 
on-shell partons (quarks and gluons) moving in the same direction (called equivalent-photon 
or Weizsacker- Williams approximation in QED (23J). For such a beam of gluons, the helicity 
is a good quantum number. The total beam helicity — the difference between the numbers 
of particles with positive and negative helicity — is also a physical observable. Thus, we have 
the gauge-invariant gluon helicity distribution Ag(x), 

Ag(x) = g + (x)-g_(x) , (24) 

where g±(x) is the density of gluons with momentum x and helicity P ■ S g — ±1. 

However, in the bound state nucleon, the gluons are virtual quanta which not only in- 
volve transverse polarization but also longitudinal polarization. They are off-shell particles 
whose wave function cannot be defined in a gauge-invariant way. Thus, the gluon helicity 
distribution formally is a matrix element of a complicated nonlocal operator [5|, 

Ay(x) = ^^y 2 J ^\PS\F + «(\n)L [Xn - fi] F + a (0)\PS) . (25) 

where F a @ = l/2e al3flu F tlu and L[x n -,o] is a light-cone gauge link. The need for gauge in- 
variance dictates that it must be the gluon field strength F^ w that appears, and the gauge 
link ensures cancellation of the gauge transformation factors U from two different spacetime 
points. In the light-cone gauge, Ag(x) becomes the matrix element of a local operator. The 
QCD expression only ensures that in the IMF and light-cone gauge, it reduces to the simple 
gluon helicity counting, 

A</(aO = J 2fc f ( 2^)2 (P\a+i(xP + , k ± )a +1 (xP + , fc ± ) - a t _ 1 (xP+, k ± )a- 1 (xP + , k ± )\P) . (26) 

In light-cone quantization with A + = 0, an on-shell gluon has two physical polarizations, 
described by the transverse components of the gauge potential A±. The matrix element is, 
however, independent of the momentum in the 3-direction and can be taken in the rest frame 
of the nucleon. However, in such a frame, the gluon parton's transverse momentum definitely 
cannot be neglected, neither can its off-shellness. For these gluons, their spin projection in 
the 3-direction is no longer gauge invariant; one cannot define the gauge-invariant gluon 
helicity in a simple fashion. Therefore, Ag(x) does not have a straightforward physical 
interpretation as the property of the bound state. 

Indeed, when integrating over all x, one gets the total parton helicity operator in a 
longitudinally polarized nucleon, 

S '™ = U^I ^ rP+F+Q (D^K-,o]^(0) , (27) 

which is a nonlocal operator, and again has no clear gauge-invariant physical interpretation. 
S g has a simple interpretation only in the light-cone gauge: it reduces to the gluon spin, 

S™ U+=o = J d 3 £(E x A) 3 . (28) 

Of course, as we have discussed in Sec. II, the same operator defined in other gauges has 
entirely different bound-state matrix elements. 
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The above discussion makes clear that gluon spin is not a gauge-invariant quantity, but 
in the light-cone gauge and IMF, the gluons in a bound state may be regarded as on- 
shell and the gluon helicity becomes a useful concept and can be measured experimentally 
(similar to electromagnetism j23[). Leaving this specific gauge and frame, this quantity 
involves longitudinally polarized gluons and no longer has a clear physical interpretation, 
though it can be defined and calculated gauge-invariantly. This motivates going beyond the 
traditional sense of gauge symmetry by introducing the gauge-invariant extension (GIE) of 
some apparently gauge-dependent quantities. A GIE of a gauge-noninvariant quantity is to 
generalize the fixed gauge result to any other gauge. There are an infinite number of GIE 
for a gauge-dependent quantity depending on the gauges in which one decides to generalize. 
The gluon spin operator in the light-cone gauge has its GIE in Eq. (127|) 



C. GIE of parton observables 

Parton physics emerges in the set-up of the IMF and light-cone gauge. In a fixed gauge, 
not all results are naturally gauge invariant, as gauge fixing often masks the gauge symmetry 
properties of both gauge- invariant as well as gauge-noninvariant quantities. This is similar 
to calculations in a fixed system of coordinates where physical quantities often do not have 
apparent properties of the underlying tensors. 

One can make GIE for any partonic operators using the same approach as in the case 
of the gluon spin, so that in any other gauges, GIE partonic operators produce the same 
matrix elements as in the light-cone gauge. 

For example, a GIE of the partial derivative in the light-cone gauge is 

id^ = iDf + J* ^-L rir]9 F +1 (f ,Ci)L [r , r] . (29) 

In this way, the partial derivative in the light-cone gauge becomes gauge invariant, and 
therefore one can talk about "gauge-invariant" canonical momentum and orbital angular 



momentum. In fact, the canonical version of the spin sum rule discussed in Refs. 24. 125 



and its scale evolution [17| can be understood as "gauge-invariant" in this sense [26l . |27 



In the same sense, one may regard all fey-dependent factorization result as gauge invariant. 



D. Problems with gauge- invariant extension 

Clearly, the GIE of an intrinsically gauge-noninvariant quantity is not naturally gauge 
invariant. As such, its usage has a number of problems: 

First, the operators resulting from GIE are in general nonlocal. Nonlocal operators do 
not usually have a clear physical meaning in general gauges, although they usually do in a 
fixed gauge. It is also difficult to perform calculations which are most practically done in 
the fixed gauge. In fact, all GIE of partonic operators cannot be calculated in lattice QCD 
because they intrinsically involve the real time which becomes imaginary after the Wick 
rotation. This, for example, is the case for the total gluon helicity AG = J dxAg(x). 

Second, GIE operators do not transform simply under Lorentz transformations. While 
local gauge-invariant operators often have simple classifications in terms of representations 
of the Lorentz group, nonlocal operators involve geometrical lines as gauge links which do 
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not have simple Lorentz transformation properties. Many gauge conditions are not Lorentz 
covariant. 

Third, there are (infinite) many nonlocal operators that have the same quantum numbers. 
For instance, one can merely change the path of a gauge link to get a new operator. All of 
these operators can mix under scale evolution. Therefore, scale-mixing calculations become 
intractable j28| . 

Finally, GIE operators are in general unmeasurable. So far, the only example is offered 
in high-energy scattering in which certain partonic GIE operators may be measured. For 
GIE operators in the Coulomb gauge, to which we will return in the next section, there is 
no known physical measurement for any of them. 

For the transverse-momentum-dependent factorization, many of the so-called factoriza- 
tion formulae are intrinsically gauge dependent. There are contributions of the same order 
that have been neglected, for instance, contributions from additional transversely polarized 
gluons. Here gauge dependence is a signal that there are legitimate contributions that have 
been left out. Thus one makes little gain by claiming the result is gauge invariant through 
GIE. 



IV. GAUGE-INVARIANT EXTENSION OF THE COULOMB GAUGE 

The Coulomb gauge is a natural gauge choice in atomic physics, which as we discussed 
in Sec. II preserves the power counting in non-relativistic effective theory. In this section, 
one considers the GIE of the Coulomb gauge, i.e., extending a Coulomb gauge result gauge- 
invariantly to other gauges. For a naturally gauge-invariant quantity, this of course does 
not generate anything new. However, if one believes that the Coulomb gauge potential A% 
is special, any operator made from it can be made gauge invariant through GIE. As we 
shall see, this formalism is just what Chen et al. have proposed as a new "gauge-invariant" 
formulation jsl, 0|. The formulation was in fact based on the well-known observation that 
in a fixed frame, any electric field E can be separated into a longitudinal part E\\ and a 
transverse part Ej_. The transverse part of the electric field and the total magnetic field 
can be generated from a transverse gauge potential A±, and the longitudinal part can be 
generated from the combinations of Au and <j), and hence the gauge freedom. The Coulomb 

gauge corresponds to choosing An = 0. Although A± is now gauge invariant, it is frame 
dependent. Therefore, gauge dependence in this formalism is traded into frame dependence. 

In this section, we will first define the GIE of the Coulomb gauge. Then we argue that 
Chen et al.'s proposal is just the same. We show that the so-called gauge-invariant anomalous 
dimension of the gluon spin operator calculated in Ref. 0] is, in fact, corresponding to the 
GIE of the gluon spin in the light-cone gauge. Finally, in the last subsection, we discuss 
why the spin and orbital angular momentum separation seems possible in optics. 



A. Coulomb gauge extension 

Let us consider how to make the GIE of the Coulomb gauge. For simplicity, we consider 
QED only and use a subscript c to denote quantities in this gauge. The gauge condition is 

V • A c = , (30) 
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and A° c becomes a nondynamical variable satisfying the constraint: 



V 2 A° = —V ■ E — —p . (31) 

Let us consider any Coulomb gauge operator O c as a function of A%. Its matrix element 
evaluated in the Coulomb gauge wave functions is (ip c \O c \4> c ). Any such operator has a GIE 
O = O c (A£) by replacing every A% in it by 

K = ^% . (32) 

which depends on the gauge- invariant field F^ u . O is now manifestly gauge invariant as A M 
only appears in the field strength tensor. If we now specialize O in any gauge by letting A 11 
satisfy a new gauge condition g, calling it O g , its matrix element is then the same as that 
in the Coulomb gauge: 

(i/> a \O a \<j> g ) = mO c \<p c ) . (33) 

Note that the states \ip g ) and \<f> g ) are correspondingly obtained in the same gauge g. In this 
sense, any operator in the Coulomb gauge O c is "gauge invariant" through O, although it 
is generally nonlocal in terms of the non-gauge- fixed potential A M . 



B. Chen et al.'s proposal 

Chen et al. claimed to have found a way to describe the spin of a gauge particle through 
a new formulation of gauge symmetry The first step in their theory is to decompose the 
gauge potential A M into two parts: 

& = A£ ure + A; hys . (34) 

The part A^ nre carries all the gauge degrees of freedom and hence has no physical conse- 
quence. The other part is "physical" because it is invariant under gauge transformations. 
In their choice, they demand A p h ys to satisfy the divergence-free condition: 

V ■ A phys = , (35) 
whereas A pum satisfies the curl-free condition: 

V x A pure = . (36) 

Under gauge transformation, one has 

^4pure ^ ^purc ^purc + V^i A p h ys ^ Apj^g Aphys • (^7) 

Therefore, one can construct gauge- invariant observables using A phys , including for example, 
the gluon spin density, E x A p h ys . 

It is easy to see that this proposal is the same as the GIE of the Coulomb gauge. If one 
is to make a calculation, the simplest way is to choose, 

A purc = V • A = , (38) 
13 



which is the same as the Coulomb gauge condition, V • A = 0. Then, A phys is the physical 
part of A^ in the Coulomb gauge which also satisfies the condition V ■ v4 p h ys = 0. The 
QED lagrangian contains only v4 p h ys , and all calculations proceed in exactly the same way 
as in the Coulomb gauge, and all results are naturally the same as those in that gauge. In 
particular, the matrix element E x A p h ys is the same as E x A calculated in the Coulomb 
gauge. 

The procedure that Chen et al. proposed is to ensure calculations in any gauge involving 
Aphy S = A — V4jV ■ A yield the same result as that in the Coulomb gauge. This is the 
significance of "gauge invariance" associated with the matrix element of operators such as 
E x A phys . In fact, the propagator (0|TAp hys Ap hys |0), considered as gauge invariant and 
physical, computed in any other gauge will yield the Coulomb gauge result: 

where rf is a temporal vector rj = (1,0,0,0). This confirms the conclusion that Chen et 
al.'s approach is the same as the GIE of the Coulomb gauge. 

Of course, as discussed above, working in a fixed frame allows one to isolate the trans- 
verse and longitudinal parts of the electric field and gauge potential, and hence allows one 
to decompose total momentum into the canonical-like contributions 0]. All calculations 
involving A± are gauge invariant. However, while the gauge part of the momentum appears 
simpler, the charged particle's momentum does not have a simple physical interpretation as 
discussed Sec. II. Furthermore, A± does not transform simply under Lorentz transformation, 
and quantities that are not naturally gauge symmetric will have complicated momentum 
dependence coming from the dynamical origin. 

C. Wakamatsu's approach 

Following the work by Chen et al., Wakamatsu [§[ proposed to generalize the procedure 
of separating the gauge and gauge-invariant parts of the potential such that one can impose 
alternative conditions on the latter and still maintain the gauge symmetry of Chen et al.'s 
decomposition. He claimed that in this way one can connect different decompositions as 
well as different frames. Unfortunately, many of the claims in the paper are incorrect. 
His procedure is not gauge invariant, and different conditions on the "gauge-invariant" 
part produce different results. Moreover, when the separation is not Lorentz covariant, the 
discussion on frame independence is also incorrect. Here we focus solely on the claim that 
the evolution of the gluon spin in this decomposition is gauge invariant His definition 
of the gluon spin is not that of Chen et al., but the GIE of the light-cone gauge. Therefore, 
his anomalous dimension result must coincide with the light-cone gauge calculation of the 
gluon spin operator, and this is exactly what he found. So his calculation does not at all 
support the claim that the gluon spin he defined is gauge invariant. 

The first step in Wakamatsu's work is also to decompose the gauge field A^ 1 into A^ me and 
Ap hys . In the paper (9|, A^ hys is defined by imposing the light-cone condition n ■ A p h ys = 0, 
where n 2 = 0. The propagator (0|T[Ap hys Ap hys ] |0) computed in any gauge just yields the 
light-cone result: 

(OITOxKCMHO) = / g^(^L) ( 9 --^±|^) . (40) 
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The gluon spin operator is constructed from ^4p hys : 

<-s P in = 2Tr[F +1 4 hys - F +2 A phys ] . (41) 

This, of course, is the same spin operator as in the light-cone gauge, which can be made 
gauge invariant by the procedure discussed in the previous section. Once this operator was 
generalized to the other gauges, Wakamatsu calculated the one-loop anomalous dimension in 
the Feynman gauge, which turned out to be the same as that obtained within the light-cone 



gauge by Hoodbhoy et al. [17|] . Nevertheless, this result is hardly novel from the GIE point 
of view. 

The GIE of the gluon spin operator from the light-cone gauge is realized by choosing 

^ hys (0 = f drr Lie-,,-] F^( v -^ ± ) L[,- r] , (42) 

where L's are light-cone gauge links. In this sense, the gluon spin operator constructed in 
Eq. (1411) is gauge invariant, and its matrix elements in any other gauge should be just the 
same as they are in the light-cone gauge. This is also true for its anomalous dimension. We 
have verified our conclusion by explicit calculations in several choices of gauges. At the lowest 
order of ^4ph ys , we calculated Ajq^ at one- loop level in general covariant and noncovariant 
gauges. All the additional integrals, with respect to the Feynman gauge result, cancelled 
each other out exactly. Therefore, Wakamatsu's result in the Feynman gauge cannot be any 
thing other than that in the light-cone gauge. 

Since Wakamatsu's calculation does not agree with Chen et al.'s result in the Coulomb 
gauge 0], he suspected that the latter was calculated incorrectly [29[. Actually, we confirmed 
the Coulomb gauge calculation. The difference is due to different GIEs. Wakamatsu's 
generalized A p ^ ys does not lead to a gauge-invariant result, which is contrary to his claim. 



D. Photon spin and orbital angular momentum in atomic physics 

It has often been claimed that the photon spin and orbital angular momenta in QED can 
be separately defined and measured, and therefore must be individually gauge invariant. 
This has served as another important motivation to look for a gauge-invariant definition of 
the gluon spin. In this subsection, we discuss the examples in optics, pointing out that this 
is possible only for optical modes with a fixed frequency. 

We first recall a bit of history about photon's angular momentum. For a circularly 



ory 

polarized plane wave, R. Beth 30j was the first to measure its spin angular momentum by 
measuring the torque exerted on the quartz wave plate it passed through. As we explained 
above, this is simply the helicity and is gauge invariant in theory. In 1992, L. Allen et al. 
pointed out that Laguerre- Gaussian laser modes also have a well-defined orbital angular 



momentum [31[. Based on this, several experiments [32H34J have been set up to observe and 
measure the orbital angular momentum of a Laguerre-Gaussian photon. 

For radiation field with e~ lu)t time-dependence, using the Maxwell equation: 

B = --V x E , (43) 
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one always has the gauge-invariant decomposition: 



J — ~ J d 3 £ [fx {E*xB + Ex B*)} 

J d^[E*((xV)E t + E*xE) . (44) 



Z£o 
CO 



The two terms may be identified as the orbital and spin angular momentum, respectively, and 
are separately gauge invariant. In particular, this is true for photon electric and magnetic 
multipoles which are often used in transitions between atomic or nuclear states. The photon 
orbital angular momentum can be defined without referring to the gauge potential at all (23| . 

It can be easily checked from the above equation that the spin equals =Fl for left-handed 
or right-handed circularly polarized light, and for linearly polarized light. In paraxial 
approximation, a Laguerre- Gaussian mode with azimuthal angular dependence of exp(i/0) 
is an eigen mode of the operator L z = —id/d<p, and carries an orbital angular momentum 
of lh. It is remarkable that the experimentalists are able to find ways to detect the effects 



of the orbital angular momentum alone [32H34 



Clearly, the above procedure only applies for a specific type of radiation field. In the 
case of QCD, the gluons in the nucleon cannot be of this type. In particular, they are 
off-shell and do not satisfy the on-shell equations of motion. The best one can do is to go 
to the IMF where the gluons appear as on-shell radiation; in this way, the gluon helicity 
and orbital angular momentum could be defined and measured "naturally" in the light-cone 
gauge. Thus we are back to the discussions made in the previous section. 



V. CONCLUSION 

This paper examines the constraints of gauge invariance on physical observables in a 
gauge theory. We advocate the standard textbook view that the naturally gauge-invariant 
(local) operators must be manifestly invariant under gauge transformations. Thus, the spin 
of gauge particles, canonical momentum, and canonical OAM are not strictly-speaking gauge 
invariant. Only in special contexts, it makes sense to talk about this. These special cases 
arise 1) in non-relativistic quantum systems such as atoms where the Coulomb gauge is a 
natural choice; 2) in radiation fields where the time variation is in harmonic oscillation (such 
as the electromagnetic multipoles and Laguerre- Gaussian laser modes); 3) in high-energy 
hadron scattering in which the gauge partons appear to have two physical polarizations and 
be on-shell. In the last case, the simple physics in the IMF and light-cone gauge cannot 
be extended to that of the bound state because the same quantity does not have a simple 
gauge-invariant interpretation in the rest frame of the nucleon. In fact, there is no gauge- 
invariant notion of the gluon spin contribution to the spin of the nucleon. The IMF frame 
result can be translated into the physics of bound state only in the light-cone gauge. 

Of course, high-energy scattering has provided great tools to learn about parton physics. 
Thus, it is possible to probe, for example, the canonical momentum and OAM of the partons 
in the IMF and light-cone gauge 26, 27||. However, they cannot be taken as gauge-invariant 



properties of bound states, contrary to many assertions made in the literature so far. Also, 
many of the TMD factorization theorems are practically model assumptions, which must be 



taken with caution [14] 



Thus, the efforts of trying to generalize the gauge symmetry from a fixed gauge calcu- 
lation is of limited use. In particular, generalizing the Coulomb gauge results as "gauge- 
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invariant" in the case of QCD appears to be neither physically attractive nor experimentally 



relevant 35, 36 
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